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1 Introduction 

Let ft be a domain is M. N , N > 2, d a locally bounded and measurable function defined in 
f2 and p a real number larger than 1. This article deals with the study of positive solutions 

of 

-div (\Du\ p ~ 2 Duj - d(x)u p ^ 1 =0 in (1.1) 

which admit an isolated singularity on the boundary of tt. It is known since the starting 
pioneering work of Serrin |11| that one of the main goals for studying the regularity of 
solutions of quasilinear equations consists in obtaining Harnack inequalities. The simplest 
form of this inequality is the following: Assume Bi r CC fi and d e L oc (B2 r ), then there 

exists a constant C — C(N,p 7 r I || d|| j^oo (b 2t .) ) ) — 1 such that any nonnegative solution u 
of 11.1 \ in Bi r satisfies 

u(x) < Cu{y) V(x,y) £ B r x B r . (1.2) 

Actually this inequality is valid for a much wider class of operators in divergence form with 
a power-type growth. Among the important consequences of this inequality are the Holder 
continuity of the weak solutions of (|1.1 I) and the two-side estimate of solutions admitting 
an isolated singularity. Among more sophisticated consequences arc the obtcntion of local 
upper estimates of solutions of the same equation near a singular point. This program has 
been carried out by Gidas and Spruck for equation 

-Am = « 5 (1.3) 

in the case N > 2 and 1 < q < (N + 2)/(N - 2) 6 , and recently by Serrin and Zhou [2] 
for equation 

-div (j-Dur -2 Duj = u q (1.4) 

in the case N > p and p — 1 < q < Np/(N — p) — 1. A third type of applications of 
Harnack inequality linked to the notion of isotropy leads to the description of positive 
isolated singularities of solutions. This was carried out by Veron ^7] for 

-Au + u q = (1.5) 
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in the case 1 < q < N/ (N — 2), and by Friedman and Veron |S] for 



di.v (\Du\ p 2 Du^j +u q = (1.6) 



when p — 1 < q < N(p — l)/(iV — p). When the singularity of u is not an internal point 
but a boundary point, the situation is more complicated and the mere inequality ljl.2 \) with 
only one function has no meaning. Boundary Harnack inequalities which deals with two 
nonnegative solutions of Ijl.l jl vanishing on a part of the boundary asserts that the two 
solutions must vanish at the same rate. For linear second order elliptic equations they are 
used for studying the properties of the harmonic measure [3] (see also £Q). For p-harmonic 
function in a ball, a sketch of construction is given by Manfredi and Weitsman |10| in order 
to obtain Fatou type results. In this article we consider singular solutions of I jl.l I) with a 
singular potential type reaction term. The first result we prove is the following: Assume dQ 
is C 2 and d is measurable, locally bounded in O \ {a} for some a £ d£l and satisfies 

\d(x)\ < C \x-a\~ p a. e. in B R (a)nfl (1.7) 

for some Cq,R > and p > 1. Then there exists a positive constant C depending also on 
N , p and Cq such that if u G C 1 (f2 \ {a}) is a nonnegative solution of [) vanishing on 
dfl \ {a}, there holds 

u{y) < u£0 < Oufo) v(a , y) e n x n s t |s| = |y| (1 8) 



Cp(y) p{x) p{y) 

where p(.) is the distance function to dtt. Another form of this estimate, usually called 
boundary Harnack inequality, asserts that if u\ and ui are two nonnegative solutions of 
Ijl.l I) vanishing on d£l \ {a}, there holds 

ui(x) u 2 (x) Cuxjx) 

< — -— < -— v{x,y) £ U x \l s.t. \x\ = \y\ . (1.9) 



Cu\{y) u 2 {y) u x {y) 

for some structural constant C > 0. Another consequence of the construction leading to 
(jl.8 j) is the existence of a power-like a priori estimate: Assume ft is a bounded C' 2 domain 
with a 6 dCl, A G is an arbitrary point and d a measurable function such that 

\d(x)\ < C \x-a\~ p a. e. in fl. (1.10) 

Then there exist two positive constants a > depending on N , p, f2 and Co, and C depending 
on the same parameters and also on A such that, any nonnegative solution u G C(f2\ {a}) n 
W^(Cl \ {a})) which vanishes on dtt \ {a} verifies 

u{x)<C - P(X 1 +1 u(A) V.TGlI\{a}. (1.11) 
\x — a\ 

The precise value of a is unknown and surely difficult to know explicitely, even in the 
simplest case when u is a p-harmonic function. In several cases the value of a is associated to 
the construction of separable p-harmonic functions called the spherical p-harmonics. Another 
striking applications of the boundary Harnack principle deals with the structure of the set 
of positive singular solutions. We prove the following: Let Q, be C 2 and bounded, a G dil 
and d satisfies \1.10 j) . Assume also that the operator v i— > — div (^\Dv\ p 2 Dvj — o?(a;)w p_1 
satisfies the comparison principle in Cl\ B e (a) for any e > 0, among nonnegative solutions 
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which vanishes on dQ \ B e (a). If u and v are two positive solutions of \) in D, which 
vanish on dfl \ {a}, there exists k > such that 

fc _1 u(a;) < v(x) < ku(x) Vx £ fl. (1.12) 

Furthermore, if we assume also either p = 2, either p > 2 and u has no critical point in Q, 
or 1 < p < 2 and d > 0, there exists k > such that 

v{x) = ku(x) Vx £ Q. (1.13) 

In the last section we give some partial results concerning the existence of singular solutions 
of equations of type Ijl.l ti and their link with separable solutions which are solution under 
the form u(x) = |a;| -7 4>{x/ \x\). If d = such specific solutions, studied by Kroll and Maz'ya 
[S], Tolksdorff Kichenassamy and Veron [7], are called spherical p-harmonics. 

Our paper is organized as follows: 1- Introduction. 2 The boundary Harnack principle. 
3 Existence of singular solutions. 4 References. 



2 The boundary Harnack principle 

In this section we consider nonnegative solutions of 

-div (\Du\ p ~ 2 Duj - d{x)u v - Y = (2.1) 

in a domain fl which may be Lipschitz continuous or C 2 . The function d, is supposed to be 
measurable and singular in the sense that it satisfies 

\d(x)\ < C \x~a\" p a.e. eO (2.2) 

for some point a £ dVl and some Cq > 0. By a solution of l|2.1 \ vanishing on dft \ {a}, we 
mean a u £ C(f2 \ {a}) such that Du £ L P {K) for every K compact, K C O \ {a} which 
verifies 

J (\Du\ p ~ 2 Du.DC + dix)^- 1 ^ dx = (2.3) 
for every £ £ C ll (f2), with compact support in f2 \ {a}. 

2.1 Estimates near the boundary in Lipschitz domains 

Let f2 be a bounded domain in ffi* with a Lipschitz continuous boundary. Then there exist 
m > 1 and ro > such that for any Q £ <9f2 there exists an isometry Iq in R w and a 
Lipschitz continuous real valued function cf> defined in R^ -1 such that 

\4>{x) -4>{y)\ < m\x~y\ V(a;,y) £ x l"" 1 , (2.4) 

and 

B2r (Q) n {x = (x',x N ) = {xi, ...,x N ) : x N > 4>{x')} = Iq (fi n B 2ro (Q)) ■ 

For any A £ i? ro (Q) D9f2, A = (a', </>(a'), r > and 7 > 0, we denote by CU,r,7 the truncated 
cone 

CU,r,7 = {a? = (x' ,xn) ■ xn > 4>( a ') 7 — a'\ < xn ~ 4>( a ')} ^ B r (A) 

The opening angle of this cone is # 7 = tan _1 (l/7) Clearly, for every 7 > m and < r < ro, 
Xq 1 (CA, r ,-y) is included into fl. Up to an orthogonal change of variable, we shall assume 
that Xq = Id. We denote also by p(x) the distance from x to dil. The next result is a 
standard geometric construction which can be found in 
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Lemma 2.1 Let Q G dfl and < r < 7'o/5 and ft, > 1 an integer. There exists Nq G N de- 
pending only on m such that for any points x and y in QnB 3r / 2 (Q) verifying min{p(a;), p(y)} > 
r/2 h , there exists a connected chain of balls B\, ■■■Bj with j < NqIi such that 



e B u yeBj, Bid B i+ i + and 2Bi c B 2r (Q) H O for 1 < i < j - 1. 



(2.5) 



Lemma 2.2 Assume u is a nonnegative solution of \) in B2 r (P) where \a — P\ > 4r. 
Then there exists a positive constant c\ depending on p, N and Cq such that 



u(x) < Cl u{y) V(z, y) G B r (P) x B r (P). 



(2.6) 



Proof. By a result of Trudinger Th. 1.1], if u is a nonnegative solution of l|2.1 \ in 

such that 



B2r{P) there exists a constant C" depending on N, p, r and |d| 1/,p 
lEU holds. Furthermore C" < C 2 exp (c 3 r 



i«(iJ 2r (F)) 

on N and p. This implies (|2IT|) since, by (|^2~|) . r \d\ 1/p 
constant depending on p and Cq . 



where C 2 and C3 depend 



L°°(B 2r (P)) 



remains bounded by a 

□ 



Up to a translation, we shall assume that the singular boundary point a is the origin of 
coordinates. 



Lemma 2.3 Assume £1 is as in Lemma Wl\ with Q ^ G dft, < r < min{r , \Q\ /4} and 
u is a nonnegative solution of \2.1 )) in £> 2r (Q). Then there exists a positive constant c 2 > 1 
depending on p, N , m and Cq such that 

u(x) < c%u(y), (2.7) 

for every x and y in B 3r /2(Q) H ^ such that mm{p(x), p(y)} > r/2 h for some /i£N. 

Proof. By Lemma T2 . 1 1 there exists AT G N* and a connected chain of j < N$h balls i?i with 
respective radii and centers 2:$, satisfying (|2.5 |) . Thus 

maxu < cirninu \fi=l,...,j, (2-8) 
by the previous lemma. Therefore 1)2.7 I) holds with c 2 = . □ 



Lemma 2.4 Let < r < \Q\ /4 and u be a nonnegative solution of ^2.1 }) in -B 2r (Q) fl £1 
which vanishes on dQ n B 2r (Q). If P £ dQ (~l B r (Q) and < s < r/(l + m) so t/iai 
B S (P) C B2r(Q), there exist two positive constants 5 and C3 depending on N, p, m and Cq 
smc/i that 

u{x) < c J X ~P M S! p(u) (2.9) 
s" 

for all x G B S (P) fl f2, w/iere M s _ P (u) = max{u(z) : 2: G S S (P) fl O}. 

Proof. Since <9f2 is Lipschitz, it is regular in the sense that there exists 9 > 0, s± > such 
that 

meas (f2 c n S s (t/)) > Omeas (B s [y)), Vy G dQ,, V0 < s < si. 
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By [151 Th. 4.2] there exists S G (0, 1) depending on p, N, Co, 9 and sx, such that for any 
y G Oil, there holds 

< C f — J V(z,z')(;B s (!/)n!lxB s (j)n!l, (2.10) 

where C depends on p, N, Cq and sup B ( y ) n!2 u = M ai)1 ,(it). Because the equation is 
homogeneous with respect to u, this local estimate is invariant if we replace u by u — 
u/M Sli y(u). Thus the dependence is homogeneous of degree 1 with respect to M auy (u), 
which implies 

/ s \ 5(1-7) 

\u{z)-u{z')\ <c f — J Af 4l , B («) v(«,«')eB,(»)nnxB,(i/)nn. (2.11) 

Taking z' = P = y, s = \x - P\, we derive (|2"1T|) . □ 

If X G i?2ro(O)n0ri and r > 0, we denote by A r (X) the point with coordinates (V, + 
r). The next result is the key point in the construction. Although it follows £Q, we give the 
proof for the sake of completeness. 

Lemma 2.5 Let < r < min{2r , \Q\ /8, Sx, 2 5 } and u be a nonnegative solution of |) 
in B 2r (Q) H fl which vanishes on dfl n B 2r (Q). Then there exists a positive constant C4 
depending only on N, p, m and Cq such that 

u{x) < c 4 u{A r/2 {Q)) Wx G B r (Q) n Q. (2.12) 

Proof. The proof is by contradiction. We first notice from (|2.9 { that if P G B2 r (Q) H <9il 
verifies P S (P) n C P 2 r(<2) H and if c 5 = {2c 3 ) 1/s , there holds 

M s/C5 ,p( w ) < ^M. iP (u). (2.13) 

By Lemma l2~31 if y G B 3r / 2 (Q) satisfies u(y) > c 2 UA r/2 (Q), it means that p(y) < r/2 h . Let 
M > such that 2 M > c 2 (defined in Lemmal2"3)l. TV = max{l+E(6+M lnc 5 / In 2), M + 5}, 
so that 2 W > 2 6 c^, and C4 = c 2 . Let u be a positive solution of 1)2.1 ll vanishing on 
B 2r (Q) (~1 <9S1 which satisfies 

u(Y ) > c%u(A r/2 (Q)), (2.14) 

for some Y Q G B r (Q) n fi. Then p(y ) < Let Q Q G <90 such that p(y ) = |*o - Qo\- 

Then 

|Q - Qo| < \Y - Qo\ + \Y -Q\< r/2 N + r < r(l + 2~ 5 ) 

Therefore Q G B 3r/2 (Q) n 90. Set s 2 = cfr/2 N , then // ;()„: C P r (i+ 2 -=+ 2 -e) (Q) C 
B 3r / 2 (Q) because s 2 < 2~ 6 by the choice of TV. Applying l|2.13 ll with s = s 2 yields to 

M S2 , Qo (u) > 2 m M S2/c m iQq (u) > 2 M u(Y a ) > 2 M c^u(A r/2 (Q)) > c^ +1 u(A r/2 (Q)), 



since \Y — Qo\ < r/2 N — s 2 /cf f and 2 M > c 2 . Hence we can choose Y\ G B S2 (Qo) n fi 
which realizes M S2 ^q (u) and this implies that p(Yx) < r/2 N+1 . A point Qi G <9f2 such that 
p(Yx) = \Yx — Qx\ satisfies also 

\Q - Qx\ <\Q- Qo\ + \Qo - Qx\ < r(i + 2- 5 + 2- 6 ). 
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Now 

M S2/2 . Ql {u) > 2 M M r/2N+ ^ Ql (u) > 2 M u(Y 1 ) > 2 2M c?u(A r/2 (Q)) > 2 N+2 A r/2 (Q). 

Iterating this procedure, we construct two sequences {Y k } of points such that p(Y k ) < 
r/2 k+N and{Q fc } such that Q k e dVl and \Q - Q k \ < r(l + 2(2~ 5 + 2~ 6 + ... + 2- 5 - k )) < 3r/2 
and 

u{Y k ) > 2 N+k A r/2 (Q) Vfc e N*. 

Since Y k e B 3r / 2 and p(Y k ) — > as k — > oo we get a contradiction with the fact that 
A r/2 (Q) > 0. □ 

Remark. The proof of the previous lemma shows that estimate l|2.12 ll is valid for a much 
more general class of equations under the form 

—divA(x, u, Du) + B(x, u, Du) = (2.15) 

where A and B are respectively vector and real valued Caratheodory functions defined on 
il x R x M. N and verifying, for some constants 7 > and eto, 01, Co > 0, 

A(x,r, q).q > j\q\ P , 

\A(x,r,q)\ ^aokf^+a.lrr 1 , 

and 

\B(x,r,q)\ <C |rr 1 |x|- p ) 

for {x,r,q) 6 x M x M w . 

2.2 Estimates near the boundary in C 2 domains 

From now we assume that fHs a bounded domain with a C 2 boundary. For any x € <9f2, 
we denote by v x the normal unit outward vector to d£l at x. Let Rq > be such that 
for any x € dCl, the two balls Br (x — Rqv x ) and Br (x + R$v x ) are subsets of and f2 c 
respectively. If P 6 dtt, we denote by N r (P) and N r (P) the points P — r^p and = P + r^p. 
Notice that r < R implies p(N r (P)) = p(Af r (P)) = r. 

Lemma 2.6 Let Q £ dfl \ {0}, < r < min{Po/2, |Q| /2} and u be a nonnegative solution 

<>f ED 

m B 2r (Q) (~l fi which vanishes on B 2r (Q) n <9f2. T/ien £/iere exist 6 e (0,2/3) and C6 > 
depending respectively on N , p and Co and N , p, Rq and Co such that 

t u(N t (P)) c 6 t . , . 

< , , < — (2.16) 

c 6 r u{N r / 2 (Q)) r 

for any P € P r (Q) n dfi and < t < rb/2. 

Proof. Up to a dilation, we can assume that \Q\ = 1, since if we replace a: by x/ \Q\, equation 
(|2.1 I) and the estimates <|2.16 I) are structuraly invariant (which means that Co and the C 
are unchanged), while the curvature constant Po is replaced by Rq/ \Q\ which is no harm 
since Q is bounded. 

Step 1 The lower bound. For a > and a > to be made precise later on, let us define 



e -a(\x-N r/2 {P)\/rr _ e -a/2° 
v(x) = V(\ X - N r/2 (P)\) = - 
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for x G B r/2 (N r/2 (P)) n B r/i (P). We set s = \x - N r/2 (P)\. Since \Q\ = 1, the function d 
satisfies — (% < < Co- Next 

-d^(|£tf ~ 2 Du) + Cov^ 1 = - \V'r 2 (0 - 1)V" + (iV - l)V'/s) + C Q V p - x . 
Therefore this last expression will be nonpositive if and only if 

(p - i) (2Z!f + i - a ) + l - N > Co (£) ^ eMW/ ( e -(./r)- _ e-/^)^ 1 

(2.17) 

where 8 = p + (1 — p)a. But 6* = by choosing a = p/{p — 1), thus Q2.17 I) is equivalent to 



<r-i)(^ + i-<0+i-"*&Q ('-« 

If r/4 < s < r/2 < 1/4 



while 



. / aas \ „ tact. \ pa 

(p-\)[— + l-a)+l-N>(p-l)(- + l-a)+l-N=^-N, 



Therefore, if we choose a such that 

-iV)>C , (2.18) 

we derive 

-div(\Dv\ p ~ 2 Dv) + Cqv"- 1 < (2.19) 

mB r/2 {N r/2 {P))nB r/i {P). Furthermore p(x) > r/16foranyx G dB r/4 (P)DB r/2 (N r/2 {P)), 
therefore 

u(x) > c 2 4 u(N r /2(P)Hx), (2.20) 

by Lemma [2.31 and since v < 1. Because u is a supersolution for i|2.19 () . we obtain that 
(|2~2TT|) holds for any x G B r/4 (P) n B r/2 (N r/2 (P)). Finally 

if x = iVt(P) with < t < r/2. This gives the left-hand side of (|2.16 \ . 

Step 2 The upper bound. Let b G (0, 2/3] be a parameter to be made precise later on. By 
the exterior sphere condition, B 3 br(J^3rb(P)) C fi c . Let fa be the first eigenfunction of the 
p-Laplace operator in B% \ B\ with Dirichlet boundary conditions and Ai the corresponding 
eigenvalue. It is well known that fa is radial. We normalize fa by fa(y) = 1 on {y : \y\ = 2} 
(notice that fa is radial) and set 



fa b (x) = fa 

thus 



\x-M rb (P)\ 

rb 



-div (\Dfa b \ p 2 Dfa b ) = j^firb 1 
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in £?3rb(.A/w,(P)) \ B rb (Af rb (P)) and vanishes on the boundary of this domain. For b small 
enough Xi/(rb)P > 1 + C for any r G (0, 1/2], thus 



-div (\Dq> rb \ p - 2 D4> rb ) - Co^ 1 > ^ rb ~ l 
in O n B 3rb {Af rb (P)) \ B rh (Nrb(P)) 3 O n B 2rb (J\f rb (P)) while M verifies 

-dro f |L>u| p ~ 2 Dm) - Cqv?- 1 < 



(2.21) 



(2.22) 

in the same domain. We can also take b > such that B 2br (Af br (P)) C B r (Q), thus 

u(x) < c 4 u(A^ r/2 ((g)) 

for x G dB 2rb (Af rb (P)) Oil by Lemma l2~Hl Now the function </> r b = C4u(N r / 2 (Q))<fi rb satisfies 
(l2~2T1 in ft n B 2rb (Af rb (P)) and dominates u on 5(0 n B 2rb (M rb (P))) = (dB 2rb {M rb (P)) H 
0) U (B 2rb (Af rb (P)) n 90). By the Diaz-Saa inequality gj 



nnB2rt(A/V 6 (P)) 



(i-c^r 



D<f>, 



rb 



V 



Kb 1 



K - 4> p rb )+dx < 0, 



valid because (u p - <^ 6 )+ G M^ 1,p (O n B 2rb (Af rb (P))) . Therefore 



/ 



K - f rb ) + dx < o, 



/nnB 2W ,(A^ ri ,(P)) 

from which follows the inequality u < 4>r b in SI n B 2rb (N rb (P)) ■ In particular 

\N t {P) - N rb {P)\ 



u(N t (P)) < c 4 0i 



u{N r/2 {Q)). 



Since </>i(s) < C(s — 1) for s G [1, 2], we obtain the right-hand side of (|2.16 \ . 
The main result of this section is the following 



□ 



Theorem 2.7 There exists two constants a > and c-j > depending on N, p, Co and TV, 
p, Co and Rq respectively such that if u is any nonnegative solution of 12.1 \) vanishing on 
dVl \ {0} there holds 



1 

C7 ' 



p(x) \x\ a 1 u(A) < u(x) < C7p(x) \x\ a 1 u(A) 



(2.23) 



for any where A is a fixed point in O such that p(A) > Rq. 



Proof. Step 1: Tangential estimate. Let such that |x| = 2r < Rq and p(x) = t < br/2. 

Let Q G dfl \ {0} such that \Q\ = \x\ and x G B r {Q), the previous lemma implies 



2 2cr 
" p(x)u(N r/2 (Q)) < u(x) < ^p( x )u{N r/2 {Q)). 



c 6 x 



(2.24) 



There exists a fixed integer k > 2 such that we can connect two points lying on <9i?2r (0) PI <90 
by fe connected balls Bj (j = 1, k) with radius r/4 and center on dB 2r (0). In particular 
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we can connect N r / 2 (Q) with A^r(O) = — 2tvq and all the balls can be taken such that the 
distance of their center to dfl be larger that r/2. Since by Lemma \2. 21 there holds 



sup it < ci inf it Vj = 1, ...fe, 

we derive 

- T l—p( x )u(N 2r (0)) < u(x) < ^ p ( x )u(N 2r (0)). (2.25) 
c i c 6 Fl Fl 

Let A = -R Q v , h = -2rv Q = N 2r {0), for I > 2, b e = -2(1 + 3(2 e ~ 1 - l)/2)rv and 
ri = 2 i ~ 1 r. Applying again Lemma 12.21 in B 2riL {bi) C fl, we have 

sup u < ci inf u W = 1,2,... (2.26) 

Let r be the solution of 

2(1 + 3(2 r - i - l)/2)r = i? <f=> t 



t i i , /.)>.. r» _ _ ln (^o +r) - In 3r 



In 2 

and = IE( r ) + 1, then Ao € -B rf() (be ), and the combination of l|2.25 I) and (|2.26 I) (applied 
£q times) yields to 

p{x)u(A ) < u(x) < 1 - p{x)u{Ao). (2.27) 



cf +<0 C6 \x\ F, 
Since r < Ro/2, the computation of r yields to 

2T = 2(R +r) <R»^ cl <(R^ hlCl/ln2 
3r r 

This implies 1|2.23 (I with a = lnci/ln2. 

5iep i?: Internal estimate. If a; G satisfies \x\ < i?o and p(x) > 6/4 absx, we can directly 
procede without using Lemma 1231 Using internal Harnack inequality (|2.6 I) and connecting 
x to iV r (0) and then to Aq we obtain 

— \x\ a u{A) < u(x) < c 7 \x\~ a u(A), (2.28) 

C7 

from which l|2.24 I) is derived since p(x) > b \x\ /4. Finally if \x\ > Ro and p(x) < Ro, we 
can replace the singular point by a regular point B £ dfl such that \x — B\ < Rq. The 
previous procedure leads to the same estimate. At end, if p(x) > Rq we apply again the 
internal Harnack inequality (|2.6 |) . Since fl is bounded, x and Aq can be joined by at most 
d = 2 diam (Q)/Rq balls Bi with radius Ro/2 and center bi satisfying p(bi) > Ro. Then 
using d times (|2.6 ll yields to H2.24 p . □ 



Remark. If p = 2 and d is regular, it is proved that Lemma 12.61 holds even if dfl is Lipschitz 
continuous. The previous proof is adapted by replacing the doubling property of the radius 
on the connecting balls B re (bi) by radii such that r<> + i = /3r^, where /3 > depends on the 
opening of the standard cone C associated to the inside cone property of fl. This observa- 
tion shows that in the general case p ^ 2 and d singular, the validity of Lemma 12 . 61 implies 
Theorem 12 . 71 when fl is a bounded domain satisfying the inside cone property. 
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Remark. In the case p = 2 and lim^^o M 2 d(x) = the value of a is known and equal to 
N — 1. When p =/= 2 the value of a is unknown, even in the case where d = 0. 

The next result is a consequence of the method used in the proof of Theorem 12. 71 

Theorem 2.8 Let u € C 1 ^ \ {0}) be a positive solutions of |) vanishing on B 2 r (~1 
(<9f2 \ {0}). Then there exists a constant eg > depending on p, N , Co and Ro such that 

1 u(y) u(x) u(y) ,„ „„. 

< -7-4 < C 9-TT' ( 2 - 29 ) 



eg p(y) 
for every x and y in Br q (0) n f2 satisfying \y \ /2 < |a;| < 2 |y|. 
Proof. By (|2.25 I) we have 

for any x & Q such that |a;| < Ro/2 and < 6 \x\ /4. If we assume that x £ SI fl -Br /2(0) 
verifies |x| < i?o/2 and > 6|ir /4 we can connect x to N\ x \(0) by a fixed number n of 
balls of radius b \ x\ /8 with their center at a distance to d£l larger than b\x\. The classical 
Harnack inequality yields to 

4^m(°)) < W W < c Wm(°))- 

c i 

Since p(a;) < |a;| < p(x)/b, we obtain, for any x € B Ro / 2 {0) H fi, 

-L-p^Ju^^O)) < u(ar) < ^p(oj)«(JV| x |(0)). (2.30) 
eg | £c | | a; I 

where eg depends on p, N, Co and R . By Harnack inequality, we can replace u(Ni x \ (0)) by 
u(JV a (0)) for any |x| /2 < s < 2 |x| and get 

r p(aOu(JV a (0)) < u{x) < ^-p(x)u(N s (0)). (2.31) 



ClCg \x\ 

If y S i3fl / 2 (0)nfl satisfies |x| /2 < |y| < |sg|, we apply twice l|2.31 \ and we get (|2.30 I) with 
c 9 = cfcg. □ 

Another consequence of this method and of Lemma \2. 21 and Lemma \2. 51 is the 

Theorem 2.9 There exists a constant c' 9 depending on N , p, Co and Ro such that any 
u G C 1 (Cl\ {0}) be a positive solutions of \2.1 \ vanishing on B 2Rq ^ (P^ \ {0}) verifies 

u{x) < c 9 u{N r (0)) (2.32) 

for every < r < R /2 and any iGilfl B 2r (0) \ B r/2 (0). 

Remark. Since Lemma 12.21 and Lemma 12.51 are valid in Lipschitz continuous domains and 
the construction of connected chain of balls too by Lemma l2.ll the above inequality remains 
valid if f2 is Lipschitz continuous. 

The next result is known as the boundary Harnack inequality. 
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Theorem 2.10 Let Q € dtt, < r < min{i?o/2, \Q\ /2}, and u\ and u 2 be two nonnegative 
solutions of \ 6 £.l J) in B 2r (Q) H f2 which vanish on B 2r (Q) H c?S7. T/iera i/iere exists Cio > 
depending respectively on N , p and Co such that 

1 mi (a;) u 2 (x) , ui(x) 
cwui(y) u 2 {y) ui{y) 

for any x, y G B r (Q) n fi. 

Proof. If a; £ P r (Q)nJ7 satisfies jo(x) < 6r/2, we denote by P^ = P the unique projection 
of a; on dil and put £ = p(x). By p. 16 |> - 

< HiW < C A (2 34) 

c 6 r " Ui(N r/2 (Q)) ~ r 

1 Ul (a;) < M*) < c 2 "ifr) 



c§ui(iV r/2 (Q)) " u 2 (AT r/2 (Q)) " b ui(7V r/2 (Q)) 

from which i|2.33 I) is derived with a first constant cio = c\. Next, if x € B r (Q) D $1 satisfies 
p(x) > br/2, we denote /3 = 2 + E(- In 6/ In 2), thus p{x) > r/2 13 . By Lemma IO 



1 Ui(x) 

4 " Ui(N r/2 (Q)) 



< < cf , (2.36) 



for i = 1,2. Therefore (|2.35 I) holds with c\ replaced by cf. Finally (|2.33 I) is verified with 
cio = max{c|, c 2 ,' 3 }. □ 

The next result is another form of the boundary Harnack inequality 

Theorem 2.11 Let ui £ C 1 (J7 \ {0}) (i = 1,2) be two nonnegative solutions of \2.1 )) 
vanishing on B 2 r ^ \ {0})- Then there exists en > depending respectively on N , p 
and Co such that for any r < Rq 



sup ( : x e n n (P r (0) \ B r/a (0) J < en inf f : x e fi n ( B r(°) \ B rM0) 

u 2 (x) J \u 2 (x) 



(2.37) 



Proof. Applying twice Theorem 12. 81 we get 



1 ui{x) < u 2 (x) < c2 ui(x) , 2 3g s 

eg ui(y) - u 2 (y) ~ 9 wi(y)' 

for any x and y such that |x| /2 < |y| < 2 |x|. Equivalently 

1 m(x) 2 mi(x) 
c%u 2 (x) u 2 (y) u 2 (x) 

which the claim with en = c§. □ 



2.3 The set of singular solutions 

We still assume that fl is a bounded domain with a C 2 boundary containing the singular 
point 0. We introduce the following assumption on the function d. 
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Definition 2.12 A measurable function d satisfying \2.2 ]) with a = £ dfl is said to 
satisfy the local comparison principle in Cl if, for any e > and any Ui £ C 1 (f2 e ) (i = 1,2) 
nonnegative solutions of \2.1 \) in Q e — Q \ B e (0) which vanish on d*Q e = dQ \ B e (0), 
u\{x) > U2(x) on Q n dB e (0) implies u\ > ui in Cl e . 

Clearly, if d is nonpositive it satisfies the local comparison principle. However there are 
many other cases, depending either on the value of Co or the rate of blow-up of d near 
which insure this principle. 

Theorem 2.13 Assume d satisfies the local comparison principle and there exists a non- 
negative nonzero solution u to \2.l |) in ft which vanishes on dVl \ {0}. If v is any other 
nonnegative solution of \2.1 \) in Q vanishing on dfl\{0} there exists k > such that v < ku. 

Proof. Since any nontrivial nonnegative solution is positive by Harnack inequalities we can 
assume that both u and v are positive in VI. We denote by TL the set of h > such that 
v < hu in Cl and we assume that TL is empty othcrwhile the results is proved. Then for any 
n £ N* there exists x n £ f2 such that v(x n ) > nu{x n ). We can assume that x n — > £ for some 
£ £ f2. Clearly £ £ fi is impossible. Let us assume first that £ £ dil \ {0} and denote by 
the projection of x n onto dVl. Thus 

p{x n ) ~ P{x n ) 

Because u and v are C 1 in ft \ {0}, 

v(x n ) - v(£ n ) dv u(x n ) - u(£„) du 
hm — = — (£) and hm — = — (£). 

n-+oo p{X n ) OV^ n— >oo p(X n ) OV^ 

Since Hopf boundary lemma is valid (see |14j). the two above normal derivative at £ are 
negative, which leads to a contradiction. Thus we are left with the 0. Set 

r n = \x n \. By Theorem EH] 

• f f V ( X ) I, 1 . -l V{Xn) . _1 

[u(x) J u(x n ) 

By the local comparison principle assumption, v > nc^w in f2 r „. This again leads to a 
contradiction. □ 

The next statement is useful to characterize unbounded solutions 

Proposition 2.14 Assume u is a nonnegative solution of \2.1 )) vanishing on dd \ {0}, 
unbounded and without extremal points near 0. Then 

\x\ u(x) ,„ 
lim 1 ' \ = oo. 2.40 

x^0 p[Xj 

Proof. Assume that (|2.40 I) is not true. Then there exist a sequence {s n } converging to 
and a constant M > such that 

f lacl u(x) 1 
sup | 1 ; : |x| =s„| < AT. 

Therefore sup : |x| = s n } < M. Because u has no extremal points near 0, say in B So (0) 

for some so > 0, the maximum of u in fl n (i? So (0) \ B Sri (0)) is achieved either on |x| = so 
or on |i| = s n . Therefore 

max{w(x) : x £ Q n (B SQ (0) \ B Sn (0))} < max {Af, max{u(a;) : x £ O n 9B So (0)}} = M. 
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Since this is valid for any n, it implies that u is bounded in Q, contradiction. □ 

Such a solution is called a singular solution. The next result, which extends a previous 
result in [2], made more precise the statement of Theorem 12. 131 

Theorem 2.15 Assume d satisfies the local comparison principle and there exists a positive 
singular solution u to 12.1 \) in Q vanishing on dfl \ {0}. Assume also either 1 < p < 2 and 
d>0, or p > 2, u admits no critical point in and 

l iminf M^M >0 . (2.41) 

x^Q u(x) 

If v is any other positive solution of \ s 2,.l )) in Q vanishing on dfl \{0} there exist k > such 
that v = ku. 

Proof. Let us assume that v is not zero. By Theorem 12.131 there exists a minimal k > 
such that v < ku. As in the proof of Theorem 12 . 1 31 the following holds: 

(i) either the graphs of v and ku are tangent at some £ e tt. If we set w — ku — v, then 

-Cw -Dw = (2.42) 

where £ is a linear elliptic operator and D = d(x)(k p ^ 1 u p ^ 1 — v p ~ 1 )/w. Since fcu(£) = 
> 0, D is locally bounded near £. If p > 2 and u admits no critical point in f2, C 
is uniformly elliptic [2] for details in a similar situation). Thus the strong maximum 
principles holds and w is locally zero. Since il is connected w = in fi. If 1 < p < 2, the 
strong maximum principle holds to and we have the same conclusion. 

(ii) either the graphs of v and ku are not tangent inside Q, but tangent on <9f2\{0}. Since the 
normal derivatives of ku and v at £ coincide, C is uniformly elliptic. If p > 2 the coefficient 
D is locally bounded. If 1 < p < 2 this is not the case but D remains nonnegative. In 
both case Hopf maximum principle applies and yields to dw/dv^(^) < 0. This is again a 
contradiction. 

(hi) or v < ku in f2, dv/dv > kdu/dv on dtt \ {0} and there exists a sequence {x n } C ft 
converging to such that 

lim V -^- = k. 

n-^oo u(X n ) 

Furthermore we can assume that 

v(x n ) ( v(x) 

— sup t 



u(x n ) [u{x) 

Put a n — max{it(a;) : = r n }. By Theorem 12.91 there exists c' g > depending on N, p, Co 
and Rq such that 

u(N rn (Q)) < o„ < c' 9 u(N rn {0)), (2.43) 

which implies 

max{fi(i) : r„/2 < < 2r„} < Cga„. (2.44) 

We set u n {x) = u(r n x)/a n , v n (x) — v(r n x)/a n and d n (x) = r^(i(r„x). Then both u n and 
v n are solutions of 

-div(\Df\ p - 2 Df) -d n f*- 1 =0 



in f2„ = £!/r„ and vanish on dil n \ {0}. By l|2.44 |) . u„ and u„ are uniformly bounded in 
Q n = Q n n (B 2 {0) \ B 1/2 (0)). Since dft n n (B 2 (0) \ Bi/ 2 (0)) is uniformly C 2 we deduce by 
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the degenerate elliptic equations theory [S] that, up to subsequences, u n and v n converge in 
the C} oc (Cl n n (B2(0) \ -Bi/ 2 (0)))-topology to functions U and V which satisfy 



-div 



(\Df\ p - 2 Df) -d^fP- 1 =0 



in H (~l (S2(0) \ 5 1 / 2 (0))), where H is the half space {t] g R" : r].i/ < 0} and doo is some 
weak limit of d n in the weak-star topology of L°°. Moreover, if p > 2, 1)2.41 [) . jointly with 
(12.43 I) and l|2.44 \ , implies that 

\Du n {x)\ = > 7 (2.45) 

for 2/5 < |a;| < 8/5, where 7 > 0. We put 

£ n = mf{^:|z| = |z n |=:r n }<fc 

and = x n /r n . Up to another choice of subsequence, we can also assume that £ n — > £ and 
£ n -» £ with £ = 1. Furthermore V < kU, V(Q = kU(g) and, if £ g &ff"n (_B 2 (0) \B 1/2 (0))), 

£«)<«) <c 

In this case, and more generally if the coincidence set 3 of V and kU has a nonempty 
intersection with dH n (-82(0) \ £?i/ 2 (0))), Hopf boundary lemma applies and implies that 
V = kV in the whole domain. If this is not the case we use i|2.45 I) to conclude again by the 
strong maximum principle that V = kU in H n (^(0) \ B 1 / 2 (0))). Therefore i — k and for 
any e > there exists n e g N such that n > n e implies 

(fc - e)tt(i) < v(x) < ku(x) Vx g Q. n 5B r „(0). 

By the local comparison principle the same estimate holds in f2 rn . Since this is valid for any 
n and any e, we conclude that v — ku. □ 



3 Existence of singular solutions 

3.1 Separable solutions 

The existence of N-dimensional regular separable p-harmonic functions associated to cones 
is due to Tolksdorff ^3]. Extension to singular function is proved in ^S]. These solutions 
are obtained as follows: Let (r, a) g R+ x S 1 ^^ 1 be the spherical coordinates in R N and 
S C S 1 ^^ 1 a smooth spherical domain. Then there exists two couples (7s, ips) and (@s, 4>s), 
where 7s and (is are positive real numbers and ips and <ps belong to C 2 {S) and vanish on 
dS , such that 

U s =r" ,s tps and V s = r~P s (f>s, (3.1) 

are p-harmonic functions in the cone Cs — {(r, cr) : r > 0, a g S}. These couples are unique 
up to homothety over ips and <p>s- Furthermore the following equation holds 

-div ((aV + |Vt7| 2 )^- 2 )/ 2 Vt7) = A(a)(aV + | V^ 2 )^ 2 )/ 2 ?? in S 

(3.2) 

77 = on dS. 
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where A(a) = a(a(p — 1) + p — N) if (a, rj) = ((3s, (j>s), and X(a) = a(a(p — 1) + N — p) if 
(a, 77) = (-fs,if>s)- 

If p 2 and N ^ 2, 75 and (3 S are unknown except if S = S^ 1 = S^ 1 n {x = 
(x',xn) : xn > e in which case 75 = 1 and ^5 = ir/v- If N = 2, equation 1)3.2 [1 is 

completely integrable and the values of the 7s and (3 s are known ([H], |7j)- When p ^ 2, 
the existence of solutions to 1)3.2 t is not easy since this is not a variational problem on S. 
Tolksdorff's method (H3|) is based upon a iV-dimensional shooting argument: he constructs 
the solution v of 

( -diu (j-Dup 2 Du) = in C s = C s n{x: \x\ > 1} , g ^ 

\ w= (2-|a|) + on <9C£. 

Then he proves, thanks to an equivalence principle, that the function v stabilizes at infinity 
under the asymptotic form v(x) « |a;| (f>(x/ \x\), with (3 > 0, which gives 1)3.2 [1 and the 
function Vs- The domain S characterizes the exponent (3. The same argument applies if 
(|3.3 I) is replaced by 

-div (\Dv\ p ~ 2 Dv) + ^f-p- = in C s = C s n {x : \x\ > 1} 

V ' W (3.4) 

v = (2- |ar|)+ on 0C£. 

with c > 0. This gives rise to a solution of l|3.4 \ in C5 under the form Vs, c = r - ^ c ' s fy where 
(3 c .s > and 

-di« ((/3 C 2 S ?7 2 + |V77| 2 )(p- 2 )/ 2 Vt7) + erf- 1 = A(/3 c , s )(/3 2 s 77 2 + I Vr?! 2 )^ 2 )/ 2 ?/ in 5 

77 = on dS. 

(3.5) 

With these considerations we can construct singular solutions of l|2.1 \ under a very restric- 
tive geometry for fi, where S = S 1 ^ -1 , the upper half unit sphere. 

Theorem 3.1 Assume d(x) — —c\x\ p with c > and Q is a bounded domain with a C 2 
boundary containing 0. Assume also dtt is flat in a neighborhood of and x.vq < for 
any x 6 O. T/ien i/iere exists a positive solution of \2.1 \ which vanishes on dfl \ {0} and 
satisfies 

lim Ixf"- 3 "' 1 u(x) = 77(0-) (3.6) 

x — ► 

x / I a; I — *■ ct 

uniformly for a G S*^ -1 , where rj is a positive solution of \3.5 \) . 

Proof. Since f2 is located on one side of the hyperplane A = {x : x.vq = 0}, the restriction 
to 57 of the function V : x *—>■ \x\ c S + T\(xj \x\) is a singular solution of 



-div 



(\Dv\ p - 2 Dvj + \v\ p - 2 v = (3.7) 



which vanishes on A n dft \ {0} and is positive on f2 \ H. Let K = max{V c s n-i(x) : x G 
dVl \ A}. Then V c s «-i — K is a solution (or a subsolution if c > 0). For any e > let u e be 
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the solution of 

-div ( \Du t \ p ~ 2 Du e ) + \u e \ p ~ 2 u e = in ft \ B e (0) 
V / \x\ 

u e = on<9ftn5 e c (0) (3.8) 
u € = V cS n-i on ft n dB € (0). 



Then V c s rr-i — K < u e < V c and the mapping e i— > u e is increasing. Therefore M e 

converges to u in the Cl oc (Cl \ {0}) topology and u is a solution of (|2.1 I) which vanishes on 
3ft \ {0} and satisfies □ 

Remark. If N = p the set of p-harmonic functions is invariant under the Moebius group, 
and in particular under the transformation x \— > 2"(x) = xj |a;| 2 which preserves Cg. In such 
a case (3 s n-i = 1. By using the transformation I it is possible to prove (see that there 
exist positive iV-harmonic functions in any bounded domain ft having a singularity at a 
point a of the boundary and vanishing on 5ft \ {a}. 

Remark. When p = 2 it is possible to prove the existence of a singular solution to 

— Au + d{x)u — in ft , , 

u = on 9ft \ {a} ^ ' 

where a G 9ft, for any C 2 domain ft and any d locally bounded in ft \ {a} such that 

— oo < liminf \x — a\ 2 d(x) < limsup \x — a\ 2 d(x) < N 2 /4. 

x >a x — fa 

We conjecture that such a result holds for (|2.1 I) and p ^ 2 although the precise upper limit 
as x — > a of \x — a\ p d(x). We believe that at least if N > p and limsup :I .^ a \x — a\ d(x) < 
((iV — p)/p) p , (the Hardy constant for W x ' v in R w ), such a singular solution do exist. 
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